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We use various methods of Borel integration to calculate the binding ground energies and masses 
of bb and tt quarkonia. The methods take into account the leading infrared renormalon structure of 
the (hard+)soft part of the binding energies E(s), and of the corresponding quark pole masses m q , 
where the contributions of these singularities in M(s) = 2m q + E(s) cancel. Beforehand, we carry 
out the separation of the binding energy into its (hard+)soft and ultrasoft parts. The resummation 
formalisms are applied to expansions of m q and E(s) in terms of quantities which do not involve 
renormalon ambiguity, such as MS mass m q and a„(fi). The renormalization scales /j, are different 
in calculations of m q , E(s) and E(us). The mass mt is extracted, and the binding energies E t i and 
the peak (resonance) energies -Ercs. for tt production are obtained. 



This is the version v2 as it will appear in Phys. Rev. D. The changes in comparison to the previous 
version: extended discussion between Eqs. (25) and (26); the paragraph between Eqs. (32) and 
(33) is new and explains the numerical dependence of the residue parameter on the factorization 
scale; several new references were added; acknowledgments were modified. The numerical results 
are unchanged. 

PACS numbers: 12.38.Bx,12.38.Cy, 12. 38. Aw 



I. INTRODUCTION 

There has been a significant activity in calculation of binding energies and masses of heavy quarkonia qq in recent 
years. The calculations, based on perturbative expansions .are primarily due to the knowledge of up to N 2 LO 
term (~ a 3 ) of the static quark-antiquark potential V(r) P, 0I and partial knowledge of the N 3 LO term there; the 
knowledge of the l/m q and 1/rn 2 correction terms (0 and references therein) and the ultrasoft gluon contributions to 
a corresponding effective theory N 3 LO Hamiltonian and the knowledge of the the pole mass m q up to order 

~ a s @> @ ■ Another impetus in these calculations was given by the observation of the fact that the contributions of 
the leading infrared (IR) renormalon singularities (at b = 1/2) of the pole mass m q and of the static potential V(r) 
cancel in the sum 2m q + V(r) Umill (analo gous cancellations were discovered and used in the physics of mesons 
with one heavy quar k 1121 1. Consequently, this cancellation effect must be present also in the total quarkonium mass 
M = 2m q + E qq [T3L Il4| . or more precisely, in M(s) = 2m q + E(s) where E(s) is the hard+soft part of the binding 
energy, i.e., the part which includes the contributions of relative quark-antiquark momenta |fc°|,|k| > m q a s , i.e., 
soft/potential scales (predominant) and higher hard scales (smaller contributions). In addition, the binding energy 
has contribution E qq (us) from the ultrasoft momenta regime |k| ~ m q a> 2 s . The ultrasoft contribution is not 
related to the 6=1/2 renormalon singularity, since this singularity has to do with the behavior of the theory in the 
region which includes the hard (~ m q ) and soft /potential (~ m q a s ) scales. 

In this work, we numerically calculate the binding ground energies E qq (separately the s and the us parts) and 
the mass (2m q + E qq ) of the heavy qq system, by taking into account the leading IR renormalon structure of m q 
and E q q{s), in the spirit of the works of Refs. \13L Il4j|. We combine some features of these two references: (a) the 
mass that we use in the perturbation expansions is a renormalon-free mass ^2 13 EH EJ 113 which we choose 
to be the MS mass rn q = m q (n = fn q ); (b) Borel integrations 1 1 - ll ] are used to perform resummations. However, 
before resummations we perform separation of the soft /potential (s) and ultrasoft (us) part of the binding energies, 
and apply the renormalon-based Borel resummation only to the s part. The renormalization scales used in the Borel 
resummations are /i/j ~ m q (hard scale) for 2m q , and m q a s < /i s < m q for E qq (s). The term corresponding to E qq (us) 
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is evaluated at [i ua ~ whenever perturbatively possible. Further, the Borel resummations are performed in three 

different ways: (a) using a slightly extended version of the full bilocal expansion of the type introduced and used in 
Refs. [ljjllll; (b) using a new "er-regularized" full bilocal expansion introduced in the present work; (c) using the form 
of the Borel transform where the leading IR renormalon structure is a common factor of the transform [19l IU (we 
call it i?-method). The Borel i nteg rations for both m q and E qq (s) are performed by the same prescription (generalized 
principal value PV [H HJ, |2j if ) so as to ensure the numerical cancellation of the renormalon contributions in the 
sum 2m q + E qq (s). Furthermore, we demonstrate numerically that in the latter sum the residues at the renormalons 
are really consistent with the renormalon cancellation when a reasonable factorization scale parameter for the s-us 
separation is used, while they become incosistent with the aforementioned cancellation when no such separation is 
used. The obtained numerical results allow us to extract the mass from the known T(IS') mass of the bb system 
and to demonstrate that the us contribution is the major source of uncertainty. We present also the numerical results 
for the ground state binding energy for the scalar and vector toponium tt. 

In Sec. [H] we recapitulate the calculation of the pole mass m q in terms of m q and a s (/ih), and summarize the 
bilocal method of Refs. 0,^3' with a slight extension in the renormalon-part of the Borel transform. In Sec. IIIII wc 
perform the separation of the binding ground energy into the soft/potential (s) and the ultrasoft (us) part, and in 
Sec. II VI we determine the s-us factorization scale parameter so that the renormalon residue reproduced from E qq (s) 
becomes consistent with the renormalon cancellation condition. In Sec. IIVI we further apply several methods of the 
Borel resummation to calculate E b j,(s) and E t t(s): the aforementioned bilocal method, the new "cr-regularized" bilocal 
method, as well as the aforementioned i?-mcthod, and always using in the expansions m q mass. We also estimate there 
the ultrasoft contributions to the binding energy. In Sec. we compare the obtained results with some of the results 
recently published in the literature and draw conclusions about the main numerical features of our resummation 
procedure. 



II. POLE MASS 



Here we redo the calculation of the pole mass m q in terms of the MS renormalon- free mass m q = m q (/j, = m Q ) 
and of a s (fi, MS), using elements of the approach of Ref. and the bilocal expansion method of Refs. 0, Hq . 
In the Borel integration, we choose the (generalized) principal value (PV) prescription fl9l I20L I2H |22||. The ratio 
S = (m q /m q — 1) has perturbation expansion in MS scheme which is at present known to order ~ a 3 (Ref. for 
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where L m = ln(^ 2 /m q ), while ft = (11 

independent coefficients with n/ = ri£ being the number of light active flavors (quarks with masses lighter than m q ) 
The natural renormalization scale here is [i — fi^ ~ m q (hard scale). 
Therefore, the Borel transform Bg(b) is known to order 
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It has renormalon sing ularities at b = 1/2, 3/2, 2, . . . , -1, -2, . . . [HEllll. The behavior of Ba near the leading 
IR renormalon singularity b = 1/2 is determined by the resulting renormalon ambiguity of m q which has to have the 
dimensions of energy and should be renormalization scale and scheme indepenent - the only such QCD scale being 
const x Aqcd 25] (cf. also 01)- This scale is proportional to the Stevenson scale A [2(| (cf. also H3)- The latter can 
be obtained in terms of the strong coupling parameter a(u; c%, C3, • ■ ■) = a s (u] c 2 , C3, . . .)/tt, where Cj = ft /ft (j > 2) 
are the parameters characterizing the renormalization scheme, by solving the renormalization group equation (RGE) 
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TABLE I: The MS RGE coefficients Ck = flu/Po (k = 1,2,3,4) and renormalon coefficients v and Cj (j = 1,2,3) for the bb 
(rif = 4) and ti {nf — 5) system. 
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where v — ci/(2/3o) = /3i/(2/3 2 ); the coefficients Cj (j > 2) will be taken here in MS scheme. Expansion of expression 
(0) in powers of a(/i) then gives 
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On the other hand, for the uncertainty in m q from the 6=1/2 renormalon singularity to be proportional to the 
quantity ©, this implies that the singular part of the Borel transform B$(b) around 6 = 1/2 must have the form 1 
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and B^ n '\b; /i) is analytic on the disk \b\ < 1. The MS coefficients c 2 and C3 are already known |2^,[2{j, but for c 4 
we have only estimates obtained by Pade-related methods. R ef. I3(M gives c 4 w 97. (n/ = 4); 86.(71/ = 5), and 

Ref. [U gives c 4 w 40. (n/ = 4); 70. (n/ = 5). However, the estimate of |30j is obtained from a polynomial in nj with 
estimated coefficients, where large cancellations occur between varios terms. Therefore, we will take as the central 
values the estimates of [3l], with the edges of the (±) uncertainties covering the values of [3(| 



c 4 = 40 ±60 
c 4 = 70 ±20 



(n/ = 4) 
(«/ = 5) 



(9a) 
(9b) 



Thus, c 3 can be obtained via Eqs. |ZbJ), (jHEJ: c 3 = 0.01 ±0.17 (n/ = 4); -0.20 ±0.08 (n/=5). The values of c fc 's and 
Cfe's are given in Table H| Now, the (full) bilocal method Q consists of taking in the expansion l|8a|l for the analytic 

part Bo" 1 '' a polynomial in powers of 6, so that the expansion of Bs around 6 = agrees with expansion @. For that, 
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the residue parameter in Eq. i|8a[l has to be determined. Using the idea of Refs. |32j it was estimated with a high 
precision in Refs. [Hill 113: 



where, according to (|8a|) 



777 1 

N m = ^-R s (b= 1/2) , (10) 

/Lt 7T 

i?s(&;/i) = (l-2&) 1+ ^ s (6 ;A1 ) . (11) 

In this work, in applications of the bilocal and related methods, we will use the value of N m as estimated in Rcf. 33], 
which used for Rs(b) truncated perturbation series (TPS) and Pade approximation [1/1]: 

N m {n f = A) = 0.555±0.020, (12a) 
N m (n f = 5) = 0.533±0.020. (12b) 

The bilocal expansion l |8a| l has then for the analytic part the polynomial 

, (m) u( m ) 

B<->(!, ;/ .) = lA m) + ^f>+ |^* 2 , (13a) 

y n— ' 

where, by convention, r§ = cq — 1. We can then take for B$ the bilocal formula, i.e., Eqs. (|8a(l and (|13f) with the 
expansion around 6 = 1/2 in the singular renormalon part truncated with the term 03(1 — 2b) 3 
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E — 
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fc=0 

Applying the (generalized) principal value (PV) prescription for the Borel integration 

1 [°° ±le ( b \ 

S(b) = - Re / d& exp --5-^ B s (b; M , (15) 

we obtain the pole mass m q in terms of the mass m q . The numerical integration is performed, using the Cauchy 
theorem, along a ray with a nonzero finite angle with respect to the b > axis, in order to avoid the vicinity of the 
pole (as explained, for example, in Refs. [20]). 

In Figs. (a), (b), we present the resulting (PV) pole masses of the b and t quarks, as function of the renormalization 
scale /Lt. The spurious /x-dependence is very weak. In addition, results of another method ("R" -method ) are presented 
in Figs. ^ (a), (b) . w ith the ^-dependence stronger in the low-/x region (fi/m q < 1). The R-method (applied in other 
contexts in Refs. |l9t |2j|) consists in the Borel integration of the function (|ll|l 

S = ^ r ±ie db exp (->) « ■ (16) 

where for Rs(b) the corresponding (NNLO) TPS is used. When we take nib = 4.23 GeV and m t = 164.00 GeV 
and we vary the values of the residue parameter N m according to Eqs. I|12l) . the bilocal method gives, at (J,/rn q = 1, 
variation Smi, — =f3 MeV and 8m t — ±20 MeV. When the central values of N m (|12|) are used, the variation of 
the obtained values of m q with fi, when (j,/m q grows from 1.0 to 1.5, is about 5 MeV and 6 MeV for nib and m*, 
respectively (for i?-method: 4 MeV and 6 MeV). When C4 is varied according to ©, the variation is about =p2 and =fl 
MeV for rrib, mt, respectively. The uncertainty in a s can be taken as a s (m T ) = 0.3254 ± 0.0125 [34|, corresponding 
to a s (Mz) — 0.1192 ± 0.0015. This uncertainty is by far the major source in the variation of the pole masses: 
(5m b ) as =t{% MeV for bilocal method (±\H MeV for i?-method), and (Sm t )a s =±m MeV for bilocal method {±1% 
MeV for R- method). 

The natural renormalization scale fi here is a hard scale /z~r?T g , and will be denoted later in this work as fi m in 
order to distinguish if from the "soft" renormalization scale \x used in the analogous renormalon-based resummations 
of the (hard+)soft binding energy E qq \s) (m q > fi > m q a s ) in Sec. IIVI The fact that the two renormalization scales 
are different does not affect the mechanism of the (6= 1/2) renormalon cancellation in the bilocal calculations of 
the meson mass (2m q + E qq (s)), because the renormalon ambiguity in each of the two terms is renormalization scale 
independent ~A, as seen by Eqs. On the other hand, if i?-type methods (fl7)f [cf. also Eq. flH| ] are applied 

for the resummations of 2m q and E qq (s), the renormalon ambiguities are renormalization scale independent in the 
approximation of the one-loop RGE running, and the renormalon cancellation is true at this one-loop level. 
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FIG. 1: The (PV) pole mass of the (a) bottom and (b) top quark, as function of the renormalization scale fi. The input parameters used 
were mj, = 4.23 GeV, rnt = 164.00 GeV, respectively; th e re sidue parameter values 1121 were used for the bilocal method. The reference 
value for a s (in MS) was taken to be a 3 (m T ) = 0.3254 (' 1341 ') corresponding to a a (Mz) = 0.1192. 



III. SEPARATION OF THE SOFT AND ULTRASOFT CONTRIBUTIONS 

The perturbation expansion of the (hard + soft + ultrasoft) binding energy E qq of the qq heavy quarkonium vector 
(5 = 1) or scalar (5 = 0) ground state (n = 1,^=0) up to the N 3 LO 0(m q a.g)was given in [33, where previous 
results of Ref. M were used. The latter reference used in part the results of Refs. 0,133,133,133 (static potential) and 
of Refs. |39l 14(1 l4lL I42I l43j (binding energy). Ref. (and 0|) employed the method of threshold expansion where 
the integrations were performed in (3 — 2s) dimensions. The reference mass scale used was the pole mass m q . The 
ground state energy expansion has the form 
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The expressions for the coefficients kij of perturbation expansion 1)17(1 for the ground state binding energy of the 
quarkonium (n = 1; I = 0; S = 1 or 0) are given below. The NLO and NNLO terms were obtained in Refs. |4oL l4ll 

their Eqs. (6) and (12), but now written in numerically more 
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Here, 03 is the hitherto unknown three-loop contribution coefficient to the QCD static potential V qq (r), whose values 
have been estimated by various methods in Refs. [l3L Il4l 133. Ei) . We will use in this work the estimates of Ref. [33| . 



obtained from the condition of renormalon cancellation in the sum (2m q 

1 o 3 ( n/ =4) « 86. ±23. , 
= 5) = 62.5 ±20. . 
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The coefficients (|19|l in the expansion (I17|) originate from quantum effects from various scale regimes of the partici- 
pating particles: (a) the hard scales (~m q ); (b) the soft and potential scales where the three momenta are |q| ~m q a s 
(\q°\ ~ m q a s in the soft and |g°| ~ m q a 2 in the potential regime); (c) ultrasoft scales where and |q| are both 
~m q a 2 . The coefficients are dominated by the soft scales; the hard scales start contributing at the NNLO and 
are numerically smaller. For this reason, in this work we will usually refer to the combined soft and hard regime 
contributions to the binding energy as simply soft (s) contribution E qq (s). Strictly speaking, it is only the pure soft 
regime that contributes to the 6=1/2 renormalon. However, for simplicity, in our renormalon-based resummations 
we will resum the hard±soft contributions E qq (s) together, not separately. This will pose no problem, since the 
hard regime, being clearly perturbative, is not expected to deteriorate the convergence properties of the series for 
Eqq(s). The natural renormalization scale \i in the resummations of E qq (s) is expected to be closer to the soft scale 
(m q a s < n < m q ). 

On the other hand, the N 3 LO coefficient A^o obtains additional contributions from the from the ultrasoft (us) 
regim e. The leading ultrasoft contribution comes from the exchange of an ultrasoft gluon in the heavy quarkonium 
[5|,[T(|. It consists of two parts: 
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where Lf « -81.538 is the QCD Bethe logarithm - see Refs. 00. 



2. The non-retarded part can be calculated as expectation value of the us effective Hamiltonian 7i us in the Coulomb 
(i.e., leading order) ground state |1), where H us (in momentum space) was derived in Refs. 0|. Direct 
calculation of the expectation value, here in coordinate space, then gives: 



— 5-^3 o(us, nonret.) 



9 



1 



47r 5 m q a 5 (/i) 



{l\H us \l) 



27tt 3 o. ,. ,1, . 

— z — a _ 1) 

8 r 



-14.196 



In 



/'■/ 



m q a 2 (fi) 



Yin 



0.9511 



7r 5 ra g a 4 (/i) 
4^ 

3 



i4 



iln 

2 (23? )» 



--ln2 
6 



^) ( l| 5(r) |l )+37 r^(l|{A r ,i}|l)} , (22a) 
mi mi r > 



(22b) 



Here, E± = —(A/9)m q a 2 s (pi) is the Coulomb energy of the state |1), and fif is the factorization energy between 
the soft (~ m q a s ) and ultrasoft (~ m q a 2 ) scale. 

In Ref. 0|, the authors included in the ultrasoft part of the Hamiltonian additional terms 5TL US which contained 
contributions from the soft regime. These terms arised because of their use of a method called threshold expansion 
|45| where the integrations over potential momenta are not performed in three dimensions but in (d — 1) = (3 — 2e) 
dimensions. However, their method gave in the soft regime also the same additional terms, but with negative sign 
(including logarithmic terms not associated with IR-divergent integrals - unphysical). Since they were interested in 
the total sum of contributions from various regimes, the method gave the correct result, as emphasized by the authors 
there. 

The s-us factorization scale /j,/ can be estimated as being roughly in the middle between the s and us energies on 
the logarithmic scale [3^1 



/'•/ 



(£ s £us) 1/2 l = k m 9 a s (/i s ) 3 / 2 , 
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where K ~ 1 and fi s w Es (< A 1 )- Therefore, the ultrasoft part of the N 3 LO coefficient fc3 ; o can be rewritten, by 
Eqs. (|21fl. (|22|l and H23|) . in terms of the s-us parameter k as 



Tk 3fi {us) = 27.512+ 7.098 ln(a,(/x a )) - 14.196 ln(«) 



(24) 
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The soft scale /z s appearing here will be fixed by the condition fi s — (A/3)rn q a s (p, s ). 

The formal perturbation expansions for the separate soft and ultrasoft parts of the ground state binding energy 
(I17JI are then 

4 2 i 3 

E qq {s) = ~-m q Tr 2 a 2 (p){l + ^2a l (fi)Y,h,jL p (fiy + a 3 ^) ^ fc 3j + a 3 (fi) [k 3 . -k 3 .o(us)} + Oia' 1 )} , (25a) 

i=l j=0 3 = 1 

E qq (us) = -^m^ 2 o 2 ( M ){a 3 (M)fc3,oM + 0(a 4 )} . (25b) 

The energy E qq (s) l)25a|) contains the leading IR renormalon effects, and E qq (us) (|25bjl does not. In these expressions, 
the common factor is the soft scale /%,(//) = (4/3)m q a s (/Lt) which is also present as the reference scale in the logarithms 
L p (/x) = ln(/i//z p (/i)) appearing with the coefficients kij (when j > 1) in Eqs. (|17|) . (JTSJ- This soft scale is equal to 
2/aB.(A*) where ae. is the (Bohr) radius of the heavy quarkonium. The renormalization scale /z in Eq. (|25a|l is of the 
order of the soft scale or above. We will re-express m q everywhere in E qq with the renormalon- free mass fn q and 
will consider the dimensionless soft-energy quantity E qq (s)/m q . 

The expansion of E og (s)/rn„ = r n (/j,)a2 +2 (^) has at large orders the seemingly peculiar feature of the 
so-called "power mismatch" (see also |4(j): when this sum is added to the expansion 2m q /m q — 

[2 + (8/3) r n {n)a n+1 {n)} , the coefficient r n (n) at powers a n (fi) of E qq {s)/m q must be combined with 
the coefficient (8/3)r„(//) at powers a n+l (ix) of 2m q /m q to ensure the cancellation of the 6 = 1/2 renormalon 
contributions. This is so because the coefficient r„(/x) contains a polynomial of n'th grade in ln[/i/(m g a(/i))] 
[cf. Eqs. H17|) . I|18|) . (|25a(l ] which, at large order n and in the large-/3n approximation, sums up approximately to 
a term ~ ((3 /2) n n\ exp\n[n/(m q a(fi))] — (f3o/2) n n\(fi/rn q )l/a(ii) [T||, effectively reducing the power of a(p.) in 
E qq (s)/m q by one. Further, the factors (/?o/2) n , n! and n in the approximate sum of the logarithmic terms in r„(/z) 
reflect the effect of the leading (6=1/2) IR renormalon in E qq (s)/m q . 

For the Borel-related resummations of E qq (s), which would account for the leading IR renormalon structure, we 
have on the basis of these facts in principle at least two possible directions to proceed. The first direction would be 
to use the Borel transform of the expansion of E qq (s) /m q = ^n(/ i ) a ™ +2 (/ i ) where the transformation a(/i) i — ^ J> is 
performed literally with respect to all a(/i)-dependence, including the one appearing in the coefficients r n . This would 
result in a Borel transform whose power expansion around the origin would include terms b k In b with t = 0, 1, 2, . . .. 

The second direction would be to divide the considered quantity by a(jl) (=>■ E qq (s)/[m q a(Jl)]), where /T, is any 
fixed soft scale, and then consider the coefficients in the expansion of this quantity in powers of a(/x) as independent 
of a(/i), e.g., by expressing them in terms of a(ji). In the obtained expansion, the coefficients now contain powers 
of logarithms ln[a(pt)] which are considered as constant (nonvariable) under the Borel transformation a([i) i— ► 6. 2 It 
is possile to see that, at large n and in the large-/3o approximation, this is equivalent to the first approach, because 
the powers of a(fi) have been decreased by one, and the coefficients are now proportional to (/3 /2)™n!/x/a(7i) where 
the factor l/a(JT) is now formally constant and does not affect the Borel transform (except as an overall constant 
factor). The equivalence is assumed to persist when we go beyond the large-/3o approximation, in the same spirit as 
the authors of Ref. lg assume their conclusions to be valid beyond large-/3o- 

We stress that in both approaches the original expansion of E qq (s) in powers of a(p) is recovered by applying the 
Borel integration according to the standard formula l|15fl term-by-term to the expansion of the Borel transform around 
6 = 0. 

In this work, we decide to follow the second direction. The main reason for this is of practical nature: The first 
approach would generate in the expansion of the Borel transform around 6 = the terms containing In 6, In 2 6, • • • , 
which introduce, at any finite order at least, a cut-singularity along the entire negative axis in the 6-plane. We are 
working at finite orders. This cut would seriously hamper our re-summations. For example, the quantity analogous to 
i?s(6) Eq. ijTTJl of the previous Section, but this time for E qq (s)/m q with the first approach, has a cut along 6 < 0, i.e., 
starting already at the origin, and the resummation at 6=1/2 would be difficult. On the other hand, the analogous 
quantity R(b) for E qq (s)/\m q a(J2)] in the second approach has no singularities at |6| < 1/2, and for 1/2 < |6| < 1 has 
only a cut without infinity along the positive axis. Such a quantity can be much more easily resummed on the basis of 
its expansion around 6 = 0. Nonetheless, the first approach presents an interesting alternative for which resummation 
techniques other than those presented here would have to be developed and/or applied. 

Thus, we will divide the soft binding energy with the quantity ~p(jl) — {A/Z)rn q a s (Ji), where Jl can be any soft scale. 
We will fix this scale by the condition Jl = (4/3)m 9 a s (/I) (=> = fi s ). Further, in the logarithms i p (/i) we express 



2 This is in close analogy wit h the behavior of t he static potential V q q(r) and its dimensionless version rV q q(r) where r ~ ae. ~ l/[m q a(^i)] 
(see, for example, Refs. [s| |To|, |Tll [T^. HI 151 ). 
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the pole mass m q in terms of m q and powers of a(/z) (cf. Sec^J, and the powers of logarithms ln fe [a(/i)] we re-express 
in terms of In [a(/i)]. This then results in the following soft binding energy quantity F(s) to be resummed 

F(s) = - A|k|>_ = a(jM) [1 + o(m)/i + a 2 ( t ,)f 2 + a 3 ( M )/ 3 + 0(a 4 )] , (26) 

where the coefficients fj depend on lna(/I) and on three scales: the renormalization scale /i (> m g a s ), the (fixed) 
soft scale ju, and m q . The coefficient /a depends, in addition, on the parameters n QEfy - fZty . p, s , and a 3 (|2U[I. The 
coefficients fj are written explicitly in Appendix 1X1 The b — 1/2 renormalon in the quantity .F (s) is then of the type 
of the renormalon of the pole mass m q discussed in the previous Sec. ITD 

However, if we divided in Eq. I|26|l by m q instead of m q and at the same time used in the resulting /^-coefficients 
lnm 9 , the numerical resummations of F(s) by methods of Sec. IIVI would give us values for E qq (s) different usually by 
not more than O(10 1 MeV) (we checked this numerically). We will briefly refer to these approaches later in this Section 
as "pole mass" approaches. A version of such pole mass bilocal approach was applied in Ref. for resummation of 
the unseparated E qq (s+us). 

The ultrasoft part il25b(l . on the other hand, has no b = 1/2 renormalon. The mass scale used there should also be 
renormalon free (jn q ). The renormalization scale fj, there should be adjusted downward to the typical us scale of the 
associated process [i i— > fi us {^m q a 2 s ) in order to come closer to a realistic estimate 3 

E qq (us) w -^TO g 7r 2 fc 3i0 (ws)a 5 (/i MS ) . (27) 



IV. EVALUATION OF THE BINDING ENERGY 

In this Section we will evaluate the soft part of the ground state energy for the vector bb [T(15)] and for the vector 
and scalar it quarkonium. In addition, we will estimate the ultrasoft part of the energy, and will extract the value of 
the mass m& from the known mass of T(IS'). 



A. Methods of resummation for the soft energy 

At first we will apply the same methods as those used in Sec. [H] However, the expansion we will use for the soft 
energy quantity F(s) l|26|l is higher by one order in a(/z) than in quantity S Eq. (JJJ of Sec.[H] In the N 3 LO coefficient 
/ 3 we have dependence on the approximately known coefficient a 3 (|2U|) , and on the s-us factorization scale parameter 

k~1 Eq. i\2'.'>\) - see Appendix lAl Eqs. I|A4(1 . It turns out that, in fa (f^), the coefficient at Iiik is larger than the 
coefficient at 03/(100 x 4 3 ). On the other hand, the coefficient at ln« in the ground state expectation value of the 
static potential is about one tenth of the corresponding coefficient in the (soft) ground binding energy 

£^(5; In k- part) w -1.93 x 10 3 (m q a 4 (fi)) In k , (28a) 
(l|Fgg(r)|l)(m re-part) w -1.95 x 10 2 (m q a 4 (^)) In k , (28b) 
£ 99 -(s;a 3 -part) = (l\V qq (r)\l)(a 3 -part) « -8.77 x 10 2 (m q a 4 (fi)) ^ ~ 3 (28c) 

Since a 3 /(100 x 4 3 ) is roughly between zero and one [cf. Eq. (J2DJ], as is also Inn, Eqs. (|28|l show that the static 
potential is more influenced by the values of a 3 than by hiK, while the situation with the (soft) binding energy is just 
reversed. More specifically: (a) the static potential is more appropriate to obtain approximate values of a 3 , as was 
done e.g. in Ref. |3^| and given in Eqs. J3UJ); (b) the soft part of the binding energy E nq (s) is more appropriate to 
obtain approximate values of the s-us factorization scale parameter k. We recall that in [33, the values of a 3 l|20|l were 
obtained by requiring that the known values of the renormalon residue parameter N m (|12f> be reproduced from the 
Borel transform of the static potential function rV qq (r). Here we will proceed analogously, and will obtain approximate 
values of k (I23|l by requiring that the residue parameter values l|12|) be reproduced from the Borel transform of the 
soft binding energy quantity F(s) of Eq. H26|) . 



3 The authors of Ref. [j^l employed a somewhat similar idea of using different evaluation methods for contributions to the spectra of 
heavy quarkonia from different regimes (short, intermediate and long-distance). A similar reasoning was employed, in the context of 
high-T QCD, in Ref. 
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As already mentioned, in contrast to the situation in Sec. [DJ the coefficients fj of the perturbation series l|26|) 
have some terms proportional to ln fc (a(/i*)) (k = 1,2, .. .) where /x» generically denotes fixed chosen scales ju, or u s 
- cf. Appendix ^ Here we will argue that these scales should be between hard and ultrasoft. These terms are 
considered constant, independent of a(/i), although they can be formally re-expressed in terms of ln fc a(/i). The terms 
of the type In a in the problem at hand are the leading terms of logarithms of ratios of various scales appearing in the 
problem (cf. Ref. Q), among them ln(Es/En) and ln(£ , us/£'s)- The typical hard, soft, and ultrasoft scales of the 
problem are, e.g., Eh = m ? , Es = (1/f), Etjs = E qq , i.e., quantities independent of the renormalization scale (/i). 4 
The /x-independent ratios of the type E$/ Eh and Eu$/E$ have expansions Ex/Ey — a(fi)[l + 0{a)\. The typical 
resummed value of this quantity can be written as a (/it*) where /i* is the typical scale of the quasiobservable Ex/Ey- 
This suggests that the lna(/z*)-terms in the coefficients of the perturbation series should really be somewhere between 
hard (En ^m g ) and ultrasoft (i?us ~ m q a l) scales. 

Similarly as in Eq. (|8a|l . we have 

+ (29) 

where the factor in front of the singular part was determined by the condition of renormalon cancellation of the sum 
2m q + E q q(s). We now define in analogy with Eq. (|llf> 

R F{s) (b;w f ) = (l-2b) 1+ "B F{s) (b;^u f ) . (30) 

Here we denoted, for clarity, explicitly the dependence on the factorization scale /i / . Expressions 1291) and (|30|l imply 

N m = 2 ^^P-R F(s) (b^^ f )\ b=1/2 . (31) 

The expansion of -Rf( s ) is exactly known up to ^b 2 , and approximately up to ~ b 3 (N 3 LO TPS), where the latter 
coefficient is dependent on k (and, more weakly, on 03). All coefficients are dependent also on the renormalization 
scale /i (> m q a s ). It turns out that the expansion i|31|) is significantly less convergent than the series l|ll[) (at 
b = 1/2). However, it is not clearly divergent, unless we take unreasonable values of k or fi. Theoretically, R F ^(b) 
should be a function with only a weak singularity (cut) at b = 1/2, and the nearest pole at b = 3/2 (i.e., the next 
renormalon pole of V qq (r) [5(| )• Thus, resummations such as Pade approximations (PA's) are expected to work better 
on R F ( s }(b) than on B F ^(b). The Pade approximation with the simplest pole structure for the N 3 LO TPS is [2/1], 
i.e., ratio of a quadratic with a linear polynomial in b. It turns out that R F ( S ) [2/1] (b) has physically acceptable pole 
structure |& po io| > 1 for most of the values of fi > m q a s and n <~ 1. Using this Pade to evaluate expression 1311) 
gives us predictions for the residue parameter N m reasonably stable under the variation of /1. On the other hand, the 
predicted value of N m depends significantly on the s-us factorization scale parameter k (1231) . 

In Fig.[2Ia) we show the dependence of N m on n, at a typical ("central") n value /i = 3 GeV, for the bb system. The 
known central value (|12a|l of N m is obtained by the R F ( S ) [2/1] (b= 1/2) expression at k « 0.59. InFig.^b) we present, 
for k — 0.59, the dependence of calculated N m on the renormalization scale /i. There, we include also the ([2/1]- 
resummed) curve for the case when no separation of the s and us parts of the energy is performed. In that case, the 
obtained values of N m are unacceptable. If the "pole mass" version is applied [mentioned in the second paragraph after 
Eq. with no separation of the s and us parts, the obtained values of the ([2/l]-resummed) curve remain above 

0.70 as well, thus unacceptable. The other values of the input parameters are chosen to have the bb "central" values: 
a 3 /4 3 = 86 1203); m b = 4.23 GeV; ju = 1.825 GeV (w fi s ) and a s (^;n/=4) = 0.3263(« a s (/z s ;n/ = 4) = 0.326) [from: 
a s (m T ; n/ = 3) = 0.3254, i.e., a s (M z ) = 0.1192 j34[]. For the RGE running, we always use four-loop MS /3-function 
(TPS) and three-loop quark threshold matching relations [5l|. with /ithresh. = 2m c , Im^. 

In Figs. |3] and 0] we present analogous results for the tt vector (5 = 1) and scalar (S = 0) bound state. The typical 
("central") values of the renormalization scale were chosen to be fj, = 55 GeV and 65 GeV, respectively. The s-us 
factorization parameter n values obtained were k = 1.16 (S = 1) and k = 1.10 (S = 0), so that i?p( s )[2/l](& = 1/2) 
would reproduce the known central value l)12bl) of the residue parameter N m . The other input parameters have 
the tt "central" values: a 3 /4 3 = 62.5 <(20b|) : m t = 164.0 GeV; p, = 31. GeV (w /Lt s ) and a s (Ji;n f = 5) = 0.1430 



4 A very similar phenomenon occurs in the perturbation expansion of the free energy of the high-temperature quark-g luon plasma, where 
the hard scale is the Matsubara frequency 2nT, and the soft scale is the Debye screening mass mj l4St |4'I . 
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FIG. 2: The residue parameter value N m as calculated from the soft part of the binding energy of the bottonium according to 
Eq. 131H . (a) as a function of the s-us factorization scale parameter k 123L at fi = 3 GeV; (b) as a function of the renormalization 
scale (j,, at k = 0.59. Further explanations given in the text. In Fig. (a), the known values $12&\ of N m are denoted as dotted 
horizontal lines. 
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FIG. 3: Same as in Fig. [5] but for the (S = 1) toponium. In Fig. (a), n - 
as dotted horizontal lines. 



55 GeV and the known values H2bl of N m are given 



[w a s (n s ; nf — 5) — 0.14]. The values of N m extracted when no separation of s and us is performed, are unacceptably 
high N m > 0.6 (also in the "pole mass" version: N m > 0.6). 

Variation N m = 0.555 ± 0.020 [bb, Eq. lfT^a|) ] implies K = 0.59 ± 0.19; variation N m = 0.533 ± 0.020 [it, Eq. lfl2b)l ] 
implies k = 1.16 1 Q (S — 1) and n — 1.10 1 g'33 (S = 0). If, on the other hand, 03 parameter is varied, according 



1.29 

to Eqs. (Dl, then for bb n = 0.59 T 0.06, and for it k = 1.16 + g;^ (S=l) and k = 1.10 + g;?? (5 = 0). Thus, the 
value of s-ms factorization scale parameter k is influenced largely by the allowed values (|12f> of the renormalon residue 
parameter, and significantly less by the allowed values CI3 (|20|) of the N 3 LO coefficient of the static qq potential. 
Therefore, we will consider the variations of N m (I12f> and of k to be related by a one-to-one relation, while the 
variations of 03 l|20|l will be considered as independent. 

In this way, we have the following values for the s-us factorization scale parameter: k (|23|) 



N m = 0.555 ± 0.020 
N m = 0.533 ± 0.020 



0.59 ±0.19 



= 1.16 



+ 0.31 
- 0.29 



1 10 + 039 
i,lu - 0.33 



(n f = 4, 5 = 1) , 
(n f = 5, 5=1), 
(n/ = 5, 5 = 0) , 



(32a) 
(32b) 
(32c) 



and thus we obtain the N 3 LO TPS Q26JI for the soft part of the ground binding energy. 
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FIG. 4: Same as in Fig. but for the scalar (S = 0) toponium. In Fig. (a), /t = 65 GeV is taken. 

We wish to add a comment on K-dependence of N m in Figs. 1 2141 Theoretically, the parameter N m should be 
independent of the s-us factorization scale /if and thus independent of the related parameter k However, among 
the /x/-dependent terms in R F ^(b; fj,; /if) of Eq. I|30l) . only the leading term is available. Due to this restrictive 
practical situation, the value of the residue parameter N m obtained by Eq. (|3 II) automatically possesses significant 
/^/-dependence (or: K-dependence), and the value of /if (or: k) is fixed by requiring that this leading order expression 
in (if reproduce the known value of N m . The value of /if obtained in this way must be physically acceptable 
if the procedure is consistent. This is analogous to the situation when a QCD observable S(Q) is known to the leading 
order ~a(/Lt) only. Equating such leading order expression SVnfQ; /x) with the known value of S(Q), a specific value of 
the renormalization scale u = /^ech is obtained such that Sm(Q; /Uech) = S(Q). This is the main idea of the effective 
charge (ECH) method |52|. If the procedure is consistent, the obtained renormalization scale n value hech should be 
of the order of the physical scale Q of the process associated with the observable: hech/Q ~ 1- The analogy with 
our case consists in the correspondence jjif /i, £us(~ E qq -) «-» Q, and [if (obtained) <-> /iech- 

Now that the value of k has been obtained, and consequently the N 3 LO TPS H2(j|) . we can perform the resummation 
of the soft part of the ground binding energy. The full bilocal method [T3.[T^ | can be performed as in Sec.lHJ Eqs. I|14|) 
and l|15p. However, now we have one term more in the TPS. Therefore 



w 27TTO 9 a(/<) (l-2b) 1+u ^ ^ fe! flS 



(33) 



where the coefficients C\. are given by Eqs. I|8b() and Q (c = 1), and the coefficients h\. in the expansion of the 
analytic part are now known up to oder k = 3 



h - N m I =jl^ w e c n T{ ;;i ± 1 -; } (. = 0,1,2,3) 

2tt m q a{n) ^ V{v + 1 - n) 



(34) 



Here, by convention, /o = 1 = Cq. Then the resummed quantity is obtained by taking the PV of the Borel integration of 
BF(s){b) of Eq. I|33(l . as in Sec.[H]for Bs(b) [Eq. (|15fl . integration along a ray]. The result would have some spurious /i- 
dependence. However, for the typical /i-scales m q > fj, > m q a s , the analytic part B^ s ^(b) of the Borel transformation 
in Eq. I|33fl turns out to have a problematic behavior in the following sense. When it is Pade- resummed as B^"'^[2/l](&), 
the obtained pole is almost always (for most /i's) unacceptably small in size: |6 po i c | < 1/2. Theoretically, B^ s ^(b) 
should have the nearest pole at b — 3/2 50]. Thus, B^? s ) appears to be too singular in the above bilocal approach, 

and the TPS and Pade evaluations of it would result in widely differing resummed values for the energy E q q(s). The 
reason for this problem appears to lie in the specific truncated form of the singular part taken in the bilocal method 
(|33|l . While the latter part describes well the behavior of the transform near b = 1/2, it influences apparently strongly 
the coefficients and thus the analytic part, so that no reliable resummation of that part (apart from TPS) can be 
done. In this context, we note that the series of terms J2k ^fc(l — 26) fc has no indication of convergence at b = 0, as 
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seen from the values of Cj in Tabled of Sec. |nj This problem can be alleviated by introducing in the renormalon part 
a "form" factor which suppresses that part away from b w 1/2, but keeps it unchanged at b « 1/2. If we choose for 
this factor a Gaussian type of function, we are led to the following set of "er-regularized" bilocal expressions for the 
Borel transform 
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N m — 



fi 



1 



2tt m q a{Jl) (l-2b) 1+ » 



x exp 



1 

"8a 



2 (1~26) 2 



E 



l + Ci(l-26) + 



C2 



C3 



A ) (1-2^ 



(35) 



The corrective terms l/(8cr 2 ) and ci/(8<7 2 ) in the coefficients of the renormalon part of Eq. (|35[) appear to ensure the 
correct known behavior of the renormalon part up to order ~ (1 — 2b)~ v+2 . The coefficients in Eq. (|35|l differ 
from hkS of the bilocal case (|34f) . and are determined by the requirement that the power expansion of expression 1351) 
reproduce the known N 3 LO TPS of the Borel transform of F s (I26[) . If a parameter increases (i.e., a > 1), formula 135|) 
is expected to gradually reduce to the bilocal formula l|33l) . If a — » 0, then the expansion of the Gaussian form function 
in 135|) would imply very large coefficients (> c~ 4 ) at the renormalon terms ~ (1 — 26)~ 1 ~ !y+fc (k = 4, 5, . . .). This is 
not expected to reflect the reality, because the results in Table [I] suggest that |cfc| < 1 for k — 4, 5, . . .. Therefore, 
we expect that the optimal choice of a would be somewhere between zero and one. Numerical analysis confirms this 
expectation. Namely, when a decreases from a — oo to about a rj 0.3-0.4, the value of the pole of the [2/1] Pade- 

resummed analytic part Bpff' (b) of Eq. (|35|l gradually turns acceptable (|6 po L 
renormalization scale /i varies in the interval [m q a s , m q ) (except close to /i 

with the TPS-evaluated and with Pade-evaluated analytic parts give for such ct's similar values of E qq -{s), indicating 
that the analytic part now manifests more clearly its non-singular behavior. When the value of a falls below 0.3, 
the analytic part starts showing erratic behavior again and the Borel resummation gives significantly differing results 
with the TPS- and the Pade-evaluated analytic parts. Further, the cr-dependence of the obtained soft energy becomes 
very strong for a < 0.3. On these grounds, the obtained optimal a turn out to be 



> 1) and rather stable when the 
m q a s ). Further, the Borel resummation 



a = 0.36 ±0.03 (n f =4, S = 1) , 
a = 0.33 ±0.03 (n f = 5, S = 0, 1) 



(36a) 
(36b) 



In Fig.lSfa) we present the (PV) Borel-resummed soft part of ground state energy for the bottonium (S — 1), as a 




FIG. 5: (a) Soft part of the ground state binding energy of bb, evaluated with the (PV) Borel-resummed expression l|85|l . as a 
function of the method parameter a. (b) Same as in (a), but for the toponium S = 1 system. Details are given in the text. 

function of the a parameter of method (|35l) . The results are given when the analytic part of Borel transform (|35() is 
either evaluated as N 3 LO TPS or as [2/1] Pade (PA). In addition, the two corresponding results (TPS, and PA) are 
given as horizontal lines when the bilocal method l|33|) is applied (a = oo). The values of the other input parameters 
have the same "central" values as in Figs. [21 and N m = 0.555 and C4 = 40. in accordance with Eqs. (|12al) and (|9a|l . In 
Fig-EJb) we present analogous results for the toponium vector (5=1) soft binding energy. The values of the input 
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parameters are the same as in Figs. |3] and 0] and in addition N m = 0.533 and C4 = 70 in accordance with Eqs. (|12b[) 
and l|9b(l . The corresponding curves for the toponium scalar (5 = 0) case are very similar to those of the 5 = 1 case. 

In addition to the methods (|33|l and Ij35(l employed up to now, which are mutually related, we want to emply as a 
cross check of our numerical results also a method unrelated to the (full) bilocal method. This will be the i?-method 
[T^. |2C|. where we resum the function R F ^(b;ii) l|30|) and then employ the (PV) Borel resummation as written in 
Eq. Ijl fijl (with Rf(s) instead of Rs there). Since we know the N 3 LO TPS of Rpr^ib), we can evaluate this function 
as TPS, or as Pade [2/1] (the Pade [1/2] is disfavored due to a more complicated and unstable pole structure). 
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FIG. 6: (a) Soft part E b j,(s) of the ground state binding energy of bb, evaluated with four different methods involving (PV) 
Borel resummation, as functions of the renormalization scale \i. Details are given in the text. In Fig. (b) the simple TPS results 
for E bb -(s) are included [Eq. I37t l. as well as the "perturbative" ultrasoft part E^?\us; [i) [Eq. (I38t l. 

The results for the soft binding energy E bb -(s) of the ground state of bottonium, as functions of the renormalization 
scale /i, are presented in Fig. E^a). The values of input parameters are taken as in Figs. El and a), and for the 
"a- regularized" method we take a = 0.36 according to Eq. I|36a|) (note that the i?-method does not need N m , C4, and 
a as input). For each of the three methods, we present two curves: when the analytic part is evaluated as TPS, or as 
Pade [2/1] (PA), where the role of the analytic part in the R- method is taken over by the function R F ^(b) itself. We 
observe from the Figure that the bilocal method (|33|l (a = oc) gives the TPS and PA results which significantly differ 
from each other. On the other hand, the "er-regularized" method (|35|> (a — 0.36) gives the TPS and PA results closer 
to each other. The methods ct-TPS, er-PA, and i?-PA give similar results in the entire presented ^(-interval. i?-TPS 
appears to fail at low /1 (w m b a s ~ 1-2 GeV). In Fig.HJb) we include, for comparison, the simple TPS evaluation of 
E bb (s), according to formula [cf. Eq. I|26|) ] 

F(s) (TPS) _ _ 9 _1„ = a(ji) rj + a(ji)h + a 2 {ll)h + a 3 (/i)/3 ] j (3?) 

47T T(X\j(X s 

where for N 2 LO TPS case we take f% = 0. In Fig. I^b) the same input parameters are used as in Fig. EJa). We see 
that the perturbation series shows strongly divergent behavior already at N 3 LO. In this Figure, we also included the 



"perturbative" ultrasoft part E^'(us;fj,) calculated according to [see Eqs. and (|25 

F^(us) = -A^J_ £ (p) (us;/i ) = k 3t0 a*M . (38) 

This quantity is highly /z-dependent. We return to the discussion of the us energy part in the next Subsection. 

In Figs.[7fa), (b), we present, in analogy with Fig.[BJa), the results for the vector and scalar toponium soft binding 
energy, respectively. The values of the input parameters are the same as in Figs. 01 El and in addition a = 0.33 
according to Eq. (|36bp . The comparative qualitative behavior of the results of various methods is similar as in the 
bottonium case, except that now i?-PA method appears to fail at low renormalization scales [i ~ m t a s ~ 30 GeV 
while i?-TPS maintains more ^(-stability there. 

In Fig. |^a) we present the results analogous to Fig. Ufa) (5 = 1 case), where we now include the results of the simple 
TPS evaluation (|37|l for ti. In Fig. JS[t>), we present the result for for the "perturbative" ultrasoft part E^' (us; fi) 
calculated according to Eq. I|38|) for the ti system (dashed /i-dependent line). We further include there the more 
realistic estimates obtained later in Subsection IIV Gl [Eq. 148all ]. 



14 



-3.12 



-3.14 - 



CD 



-3.18 - 



-3.2 



-3.22 



S=1 (a) 


K=1.16 

m, = 164.0 GeV 
a B (M z )=0.1192 






R-TPS 


-~- ><: - — _ / 


--.aPA 


R-PA 


. R-PA 




oTPS 






■PA 

a=~PA 


o 


oTPS 



20 30 40 



50 60 70 
|i[GeV] 



80 90 100 



> 

CD 

llT 



-3.16 
-3.17 



-3.2 
-3.21 
-3.22 
-3.23 
-3.24 
-3.25 
-3.26 



y— \ (b) 

- \ ' -.R-TPS 


S=0 k=i.io 

m, = 164.0 GeV 
a s (M z )=0.1192 - 


\ oTPS 






/ 


[ « ::>S ,R-PA 






:> A 


r, ..PA 





20 30 40 



50 60 70 
|i[GeV] 



80 90 100 



FIG. 7: Same as in Fig. EJa), but for the toponium system - (a) vector (S = 1), (b) scalar (5 = 0). Details are given in the 
text. 
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FIG. 8: (a) Same as in Fig.[7|^a), but now the results of the simple TPS evaluation 1371 are included, (b) The ultrasoft energy 
parts by different evaluations: (us; fi) by Eq. (1381 : E t i(us) values of Eq. 148aH as straight lines. The input parameters are 
the same as in Fig.|7Ja). 



B. Extraction of bottom mass 



We need to address now also the problem of evaluating the ultrasoft part E q q(us) of the ground state binding energy. 
The estimate of the perturbative part is given in Eq. I)27|l , where it was essential to take for the renormalization scale 
a us scale fi ~ /i ns ~ m q a 2 . 

For the bottonium case, this scale is below 1 GeV, the energy at which we cannot determine perturbatively ct s (/-*)• 
This indicates that in the bottonium the us part of the binding energy has an appreciable nonperturbative part. 
The lowest energy at which we can still determine perturbatively a s is /i « 1.5-2.0 GeV, giving a s (p) 0.30 — 0.35. 
Although this is a soft scale for bb, we will use this also as an ultrasoft scale. Then by Eq. (|27fl 

E m (us)M w -~m q TT 2 k 3 , (us)a 5 (ii us ) « (-150 ± 100) MeV . (39) 

y 

The nonperturbative contribution coming from the gluonic condensate is given by [53l | 

E bh (us)^ « m ^ 2 ^ (fn b a s {^ us )\ (a(^ us )G^Gn « (50 ± 35) MeV , (40) 
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TABLE 


II: The separate uncertainties 


Sm b (in MeV) for 


the 


extracted 


value of m,b from various 


sources: 


1.) its 


[SE bi (us 


•j(p+np) _ _ 


100 ± 106 MeV]; 2.) 


u, = 3 ± 1 


GeV; 3.) 


f^m — 


rm.fi + 5V 


4 ^ a- f-m 1 


= 0.3254 ± 0.0125 


[aJMz] 


i = 0.1192 ± 0.00151; 5.) 


N m = 


0.555 ± 0.020 


[k = 


0.59 ± 0.19]; ( 


3.) a 3 /4 3 = 


86. ± 23.; 7.) 


c 4 = 40. 


± 60.; 8. 


a = 0.36 ± 0.03; 9. m 


c ^0 (<5M T (m c /0) = 


= ±10 MeV). 


















us 






fl'ITl 






N m 




C4 


(J 


711 c 


cr-TPS 


-49 


+9 




A 

— 4 


— 13 




-3 


[ •'") 


Q 

— 5 


+4 


_5 




+49 


-13 






+ 14 




+2 


— Z 


+o 


-9 


+5 


a-PA 


-49 


+13 




-4 


-15 




-3 


+1 


-5 


+5 


-5 




+49 


-20 




+2 


+15 




+2 


-l 


+2 


-9 


+5 


R-TPS 


-50 


-4 




+4 


-8 




-9 


-3 








-5 




+50 


+45 




-40 


+10 




+11 


+3 








+5 


R-PA 


-49 


+3 




+4 


-11 




-4 


-2 








-5 




+49 


-20 




-40 


+ 12 




+4 


+2 








+5 



where we used m b = 4.2 GeV, and the value of the gluon condensate ((a s /7r)G 2 ) = 0.009 + 0.007 GeV 4 54\. Eqs. (jSHJ 
and (gUJ) give 

£^(ws) (p+np) w (-100 ± 106) MeV , (41) 

where the two uncertainties were added in quadrature. In addition, there are finite charm mass contributions which 
have been calculated in Ref. [^3 (based on the results of Refs. 0IEEIE3)- These contributions modify the values of 
m b and E bb , resulting in the contribution to the mass My{lS) — (2m b + E bb ) 

SM r (lS, m c ^ 0) « 25 ± 10 MeV . (42) 

The estimates (|41|l , and (|42|) then give a rough estimate of the us and m c ^ contributions to the bottonium mass 
5Mt(1S; us + m c ) w (—75 ± 106) MeV. The mass of the T(IS') vector bottonium ground state is well measured 
Mr (IS) = 9460 MeV with virtually no uncertainty j5^- Therefore, the pure perturbative "soft" mass is 

M r (lS; s) = 2m b + E bi (s) = 9535 + 106 MeV , (43) 

where the uncertainty ±106 MeV is the rough estimate dominated by the uncertainty of the us regime contribution. 
Our numerical results for E b i(s) in this Section and for m& presented in Sec. fallow us, by varying the input value 
of rfib, to adjust the sum 2m b ± E bb (\s) to the value given in Eq. (|43(l . For the soft binding energy we apply the 
"<7-regularized bilocal methods tr-TPS and <r-PA, and i?-TPS and i?-PA, with the aforementioned "central" input 
parameters: a s (M z ) = 0.1192; Jl = 1.825 GeV (« (j, s ), thus a s (/i,n/=4) = 0.3263 [a s (^ s ) = 0.326]; N m = 0.555; 
k = 0.59; a = 0.36; a 3 /4 3 = 86; c 4 (MS) = 40. For 2m b we apply the bilocal-TPS and R-TPS method, with 
renormalization scale n m /m b = 1, both methods giving us very similar results [cf. Fig. ^ (a)]- The bilocal-TPS 
method is applied for 2m b when cr-TPS and c-PA are applied for E bb (s); the -R-TPS is applied for 2m b when i?-TPS 
and i?-PA are applied for E bb (s) (the same combinations of methods will be applied in the next Subsec. TlV Gl to the 
study of toponium). The extracted values of m b = m b (u = rn b ) are then 

m b = 4.225 ± 0.054 GeV (cr-TPS) , (44a) 

= 4.220 ± 0.056 GeV (cr - PA) , (44b) 

= 4.243 ± 0.080 GeV (.R-TPS), (44c) 

= 4.235 ± 0.068 GeV (i?-PA) . (44d) 

The uncertainties are the combination, in quadrature, of uncertainties from various sources, shown in Table ^ for 
each of the four methods. In the case of asymetric uncertainties, the larger is taken. The largest uncertainty (±0.049 
GeV) comes from the us sector uncertainty ±0.106 GeV of Eq. (S). ^ n * ne case °f -R-TPS method, the variation of 
the soft binding energy with the variation of the renormalization scale is a competing source of uncertainty for m b 
(±0.045 GeV), and in the case of i?-TPS and i?-PA methods (where m t is resummed by -R-TPS) the uncertainty 
from the variation of the renormalization scale /x m in the 2mt-resummation is competing as well (0.040 GeV). The 
arithmetic average of the central values of Eqs. I|44|) gives us 



m b = 4.231 ± 0.068 GeV (our average) , 



(45) 
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TABLE III: The uncertainty 5E t j(us) and the separate uncertainties in SE t t(s) (in MeV) for the toponium S— 1 binding energy 
from various sources: 1.) 5E t t(us) [cf. Eqs. (IBM; 2.) M = 55 ±20 GeV; 3.) a s (M z ) = 0.1192 ± 0.0015; 4.) N m = 0.533 ± 0.020 
[k = 1.16 1 g|g]; 5.) a 3 /4 3 = 62.5 ± 20.; 6.) c 4 = 70. ± 20.; 7.) a = 0.33 ± 0.03. The pole mass m t is kept fixed at the value 
mt = 174.30 GeV. 





SE(us) 


5E(s) : n 




N m 


03 


C4 


a 


cr-TPS 


-100 


-7 


-105 


-21 


-4 


+5 


-17 




+100 


+8 


+109 


+23 


+4 


-5 


+29 


a-PA 


-100 


-10 


-105 


-21 


-3 


+3 


-16 




+100 


+8 


+109 


+22 


+3 


-2 


+26 


7?-TPS 


-100 


+3 


-104 


-21 


-8 










+100 


+8 


+109 


+26 


+8 








H-PA 


-100 


+7 


-105 


-18 


-6 










+100 


-28 


+110 


+20 


+7 









where we emphasize that the central value for the strong coupling parameter was chosen to be a s (Mz) = 0.1192. 
In Eq. H45|) . the uncertainty was chosen to be the second largest uncertainty in Eqs. (|44|) . The largest uncertainty 
±0.080 GeV of the i?-TPS method, was discarded because -R-TPS is the only one of the four methods which fails 
simultaneously at the low fi m (< mj,) and low fi (< 3 GeV) renormalization scales. 

C. Numerical results for the toponium 

For the binding energy of the toponium, the numerical results are obtained in the following way. First the value of 
the (PV) pole mass mt is fixed to the central experimental value m ( = 174.3 GeV [H^. For calculation of the binding 
energy, m t is an input parameter (but not m t ). When a s varies [a s (Mz) — 0.1192 ± 0.0015], the two methods of 
Sec. m [cf. Fig. n (b)], with the renormalization scale [i m — ra t , give WL t = 164.000+ q'IH GeV (bilocal method) and 
m t = 164.011+ ^lll GeV (i?-method), when m t = 174.3 GeV (PV value). The values of m t change by 0.020 GeV or 
less when the other parameters are varied (renormalization scale \x m \ N m and C4 for bilocal method; see Sec.[nJ), and 
such small variation in mt influences the toponium binding energy insignificantly 5 - by less than 0.001 GeV. 

We use as the central m t input value Wi t — 164.000 GeV to calculate E t i{s) with the four aforementioned Borel 
methods, using for other input parameters their "central" values used in Figs.|5fb),[3 a s {Mz) — 0.1192; Ji = 31 GeV 
(w n a ), thus a s {Jx,n f = 5) = 0.143 [a s (fx s ) = 0.14]; N m = 0.533; re = 1.16(5=1), 1.10(5=0); \i = 55 GeV (5 = 1), 65 
GeV (5 = 0); a = 0.33; a 3 /4 3 = 62.5; c 4 (MS) = 70. Then the resulting toponium soft energy is 



E t i(s) = —3.163 ± 0.116 GeV (-3.216 ± 0.120 GeV) (cr-TPS), (46a) 

= -3.158 + 0.115 GeV (-3.212 + 0.118 GeV) (er — PA) , (46b) 

= -3.154 + 0.113 GeV (-3.200 + 0.116 GeV) (i?-TPS) , (46c) 

= -3.159 + 0.115 GeV (-3.209 + 0.118 GeV) (i?-PA) , (46d) 



where the results are given for the vector (5=1) case and in parentheses for the scalar (5 = 0) toponium case. The 
uncertainties are combinations, in quadrature, of uncertainties coming from various input sources: 5a s , 5/i, Sa^, Sc^, 
SN m , and 5a. When Sa s is varied, the value Wit is varied as well, as described above, but otherwise it is kept fixed (at 
164.000 GeV). All the corresponding separate uncertainties SE t t(s) are given in Tables ITTT1 for 5=1 and IIVI for 5 = 0. 
The ultrasoft part E t t{us) is principally perturbative and can be estimated by formula l|27|l where the us coefficient 
is given by (|24|) . This part is more manageable than in the bottonium case, because the typical us energy now is still 
in the perturbative regime: fj, us ~ 10 1 GeV. We determine this energy by the condition 

^ us = K'Wita 2 s (n us ) , (47) 

where k' ~ 1. The value k' = 1 corresponds to /j, us w 7 GeV. Eqs. (|27H and 12411 then give for the value E t i(us) — 
-0.255 GeV (5 = 1) and -0.272 GeV (5 = 0). When we change to k' = 2 (fi us = 10.5 GeV), the values of E t i{us) go 



5 A variation Smt ± 10.0 MeV results in 8E t ^(s) = =p0.11 MeV, when all other input parameters are kept fixed. 
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TABLE IV: As Table ITTT1 but for 5 = 0. The input parameters are the same, except for a (= 65 ±20 GeV) and ft (= 1.10 
corresponding to N m = 0.533 ± 0.020). 





dE(us) 


oh(s) : (i 


a s 


N m 


03 


C4 


a 


<t-1Pb 


—110 


—0 


— 107 


—23 


—3 


+5 


— 17 




+110 


+8 


+ 112 


+26 


+4 


-5 


+31 


a-PA 


-110 


-8 


-107 


-23 


-2 


+4 


-16 




+ 110 


+8 


+ 112 


+25 


+3 


-2 


+28 


7?-TPS 


-110 





-106 


-25 


-7 










+110 


+13 


+111 


+30 


+8 








R-PA 


-110 


+9 


-107 


-20 


-6 










+110 


-27 


+112 


+23 


+6 









up by 0.100 and 0.110 for the 5=1,2, respectively. This we adopt as the uncertainty in the us sector. Therefore, we 
have by Eq. l|27| 

E t t (us) = -0.255 + 0.100 GeV (5 = 1), (48a) 

= -0.272 + 0.110 GeV (5 = 0), (48b) 

corresponding to /i us = 7.0 1 ^'5 GeV. When we take for the soft part E t i(s) the arithmetic average of the results of 
the four methods l|46|l . and combining it with the ultrasoft part (|48|l . we obtain 

E tl = -3.413 + 0.153 GeV (5 = 1), (49a) 
= -3.481 + 0.163 GeV (5 = 0). (49b) 

The two dominant contributions to the uncertainties in Eqs. (149(1 are the uncertainty from a s in the soft sector, and 
the uncertainty of the ultrasoft sector, as seen from Tables Hill and iTVl and Eqs. 

The results (|49|) are relevant for the future determinations of rn t from ti production near threshold. We recall that 
the determination of the pole mass rat has, due to the 6=1/2 renormalon singularity, an intrinsic ambiguity of order 
Aqce>, i.e., several hundred MeV, and cannot be determined from experiments with a higher accuracy. But the mass 
rnt could be eventually determined with accuracy of less than 100 MeV, as pointed out in Refs. ^(| where toponium 
mass was investigated using large-/?o arguments. The 5=1 toponium state is produced in e + e~ annihilation, while 
the 5 = state in unpolarized 77 collisions. The produced resonance is not exactly at the ground state mass value 
(2nit + E t i) because of the large decay width of the toponium [5JJ |6J| 

E xes . = 2m t + E t t + 8 r E 1 . cs . . (50) 

The shift in Eq. (|5T]|l is 8 r E ICS , = 100 ± 10 MeV j^E 0| and it is rather stable under the variation of all input 
parameters, including a s and m t . At this point, we should evaluate the sum (2m t + E t t) for a general input value of 
rnt (~ 164. GeV). The expected central values of (2mt + E t i) can be inferred from the central values of the binding 
energies (I49|) which were obtained with the choice m f = 164.000 GeV. We obtain the variation 

6(2m t + E t t) ~ ±2-09 5m t , (51) 

when only the input parameter rn t is varied around its central value 164.00 GeV, while all the other input parameters 
(a s , N m , /i TO , /i, 03, C4, a) are kept fixed at their corresponding central values. 6 At rh~ t — 164.000 GeV, the bilocal 
method gives m t = 174.300 GeV and the i?-method m t = 174.288 GeV. Thus, combining the average of this with 
relations l|51|l and l|49|l . we expect the approximate central values (2mt + E t i) = 345.175 GeV for 5=1 and 345.107 
GeV for 5 = 0, when rn t = 164.000 GeV. The uncertainties of (2mt + E t t) originate from the variation of all the input 
parameters except m t . Some of them are expected to be close to the uncertainties in Eqs. l|4*iH> given for the binding 
energies. However, they are not equal to these uncertainties of Eqs. (|49|l because the latter were obtained by keeping 
the pole mass fixed (m t = 174.3 GeV). Now, however, m t — 164.0 GeV is kept fixed, and variations of E t i and mt 



More precisely, Smt = ±100 MeV would correspond to 6(2mt + E t j) « ±208.8 MeV, of which S(2m t ) = ±210.1 MeV, SE tI (s) = r+L.l 
MeV, and SE tI (us) = ^0.2 MeV. 
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TABLE V: The separate uncertainties S[2m t + E t t(s + us)] (in MeV) for the toponium S = 1 mass from various sources: 
1.) /i ua = 7.0 + H GeV [cf. Eqs. @HJ]; 2.) fj, = 55 ± 20 GeV; 3.) = m b (l ± 0.5); 4.) a s (M z ) = 0.1192 ± 0.0015; 5.) 
N m = 0.533 ± 0.020 [re = 1.16 + 029]; 6.) a 3 /4 3 = 62.5 ± 20.; 7.) c 4 = 70. ± 20.; 8.) a = 0.33 ± 0.03. The input mass 
m t = 164.00 GeV is kept fixed. 
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03 
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-16 




+100 
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+3 
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+26 
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-100 


+2 


-9 
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+54 


-8 










+100 
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-95 
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i?-PA 


-100 
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-9 
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+100 
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TABLE VI: As Table|V| but for S = 0. The input parameters are the same, except for ^ (= 65 + 20 GeV) and k (= 1.10 ± 
corresponding to N m = 0.533 ± 0.020). 





us 


/" 


fJ,m 


a s 


N m 


03 
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cr-TPS 


-110 


-7 


+ 13 


+184 


+ 112 
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+ 110 
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-4 
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cr-PA 


-110 


-9 


+13 


+184 


+ 112 
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+1 


-17 




+ 110 
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-1 
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-110 


+9 


-9 


+183 


+75 


-6 










+110 


-27 


-95 


-198 


-90 


+6 









become correlated in the sum (2m* + E t i). More importantly, the variation of a s now changes E t i(s) and 2mt, and, 

to a lesser degree, E t t{us); the variation of N m changes k which in turn changes E t t(us) [Eqs. (|2"7jl and (|2"4"jl ] and, to 
a lesser degree, E t i{s) and 2mt- The explicit calculations give for S=l 

(2m t +%) = 345.181 + 0.253 GeV (cr-TPS), (52a) 

= 345.186 + 0.253 GeV (cr-PA), (52b) 

= 345.168 + 0.254 GeV (R-TPS) , (52c) 

= 345.163 + 0.256 GeV (i?-PA) , (52d) 



and for 5* = 



(2m t +E t t) = 345.119 + 0.263 GeV (cr-TPS), (53a) 

= 345.116 + 0.263 GeV (cr-PA), (53b) 

= 345.105 + 0.261 GeV (R-TPS), (53c) 

= 345.096 ± 0.263 GeV (R-PA) . (53d) 



Here, the resummation of the mass 2mt was performed by the bilocal TPS method in the first two cases [Eqs. (|52a|l . 
(l52bTl and (|53aJ) . (|53d) ]. and by the i?-TPS method in the last two cases [Eqs. (|53cl) . §2$ and (|5lc|) . (J533J)] - cf. SecUD 
In Tables Ivl and IVTI we give, for S'=l and S = 0, respectively, separate uncertainties in the mass (2mt + E t i) coming 
from the corresponding variations of the input parameters a s , N m , fj, m , fi, 03, C4, a and [A us . Adding them in 
quadrature, this gave the uncertainties in Eqs. (|52a|l - lj53d|) . We take the arithmetic average of the central values in 
Eqs. (|52a|l - H52d|) for 5 = 1, and of the central values in Eqs. I|53a|) - (|53d|l for 5 = 



(2m t +Eu) = 345.175 + 0.256 GeV (5 = 1), (54a) 
(2m f + E t t) = 345.109 ± 0.263 GeV (5 = 0). (54b) 
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Combining this with Eq. I|51|) and the aforementioned shift value S r E rcs , = 100 ± 10 MeV in Eq. I|50|) . this gives finally 

E TCS . = (345.28 ± 0.26) GeV + 2.09 (m t - 164.00 GeV) (5 = 1), (55a) 
= (345.21 ± 0.26) GeV + 2.09 (m t - 164.00 GeV) (5 = 0), (55b) 

In Tables El and EH we see that the major source of uncertainty is from the uncertainty 8a s {Mz) — ±0.0015, followed 
by the uncertainty of the ultrasoft sector scale Sfi us [cf. Eqs. I|48|l ] and in the cr-methods by the uncertainty in the 
renormalon residue parameter SN m — ±0.020 and in i?-methods by the uncertainty dfj, m in the renormalization scale 
for the resummation of 2m t . 

We could adopt in the ultrasoft regime a more conservative approach, allowing for the parameter k 1 in Eq. g3 not 
just to vary from value 1 up to value 2, but also to vary down to value 1/2. This would correspond to the variation 
of /i us from 7 GeV down to 4.27 GeV [thus increasing a s (/i us ) from 0.198 to 0.228, if keeping n/ = 5]. This would 
increase the uncertainties ±0.100 and ±0.110 GeV in Eqs. igHJl to ±0.260 and ±0.275 GeV, respectively. This would 
give in our results H55J) for the tt resonance the increased uncertainties 0.35 GeV (5=1) and ±0.36 GeV (5 = 0). 

The present experimental uncertainty in the pole mass is 5m t — 5.1 GeV [Hsj . corresponding to Suit = 4.86 GeV 
(provided we consider mt to be the Principal Value pole mass). This implies, according to results <|55|) . the present 
experimental uncertainty (c^rcsjcxp. = ±10.16 GeV, which is still very much above the uncertainties ±0.26 GeV (or: 
±0.36 when conservative approach in the us regime) coming from the uncertainties of the resummation methods and 
of the input parameters (other then m t ). 

In this work we did not include electroweak (Higgs) effects, which are significant in the case of the top quark. In 
Refs. j62,|6^| the 0(a) and 0(aa s ) corrections, respectively, to the relation between m t and m t mass were calculated. 
The size of these corrections significantly depends on the hitherto unknown mass Mh- For low Higgs masses Mh = 100- 
300 GeV, these corrections change the value of m t , for a given value of m t , by several percent. Inclusion of these 
effects would be important for a realistic extraction of rn t from the resonance energy of the tt production. 7 

V. COMPARISONS AND CONCLUSIONS 

In this Section we will compare our results with some of the results recently published in the literature. 

Our results for the mass Tab, Eqs. (|44(1 . I|45(l . Table ITU will be compared with those recently obtained by authors 
who either used pQCD expansions for the T(l S) resonance mass, or T sum rules. The only input parameter common 
to all these methods is a s . The comparison of the various methods is more reasonable if the same central input value 
of (MS) a s (Mz) is taken. Our central value was a s (m T ) — 0.3254 [=>■ a s (Mz) = 0.1192] since such [^ij, or similar 
pel Ir34| , values follow from the (nonstrange) semihadronic r decay data which are very precise [65| . On the other 
hand, the world average as of September 2002 is a s (Mz) = 0.1183 ± 0.0027 [6(j. Most of the authors during the 
last four years used central value a s (Mz) ~ 0.118. Therefore, for comparisons, we convert our results l|44|l to this 
central value of a s - more specifically, from a s (Mz) = 0.1192 ± 0.0015 to 0.1180 ± 0.0015. This can be easily done 
by inspecting in Table HP the column under a s , giving in Eqs. Q44af> - I|44d|) an increase in the central values of 11, 12, 
8 and 10 MeV, respectively. This gives the average 10 MeV higher than in Eq. (|45(l 

m b = 4.241 ± 0.068 GeV [average when: a s (M z ) = 0.1180 ± 0.0015] . (56) 

All the separate uncertainties given in Table ITT1 remain, of course, valid also in this translated result. In Table IVIll we 
give comparison of this result with others in the recent literature. All these results have the central value a s (Mz) = 
0.118. Wherever the central value of a s was different [l4l l35j. we performed the corresponding translation. There 
are two important numerical effects in our result. The first is the separate evaluation of the "perturbative" ultrasoft 
energy part at the corresponding low renormalization energy (< 2 GeV), Eqs. (|27|) and (|39|l . If we had not separated 
the ( "perturbative" ) ultrasoft from the soft part of the binding energy, the use of the common renormalization energy 
scale \i (~ 3 GeV) in the resummation then would have given us the central value of E t t{us) by about +100 MeV 
higher. Then the extracted value of nib would have gone down by about 46 MeV, giving the value nib ~ 4.195 ± 0.068, 
with the central value close to that of L03 in Tabic IVIII On the other hand, that latter value is quite clearly lower 
than the value PS02 in Table IVTTI by about 150 MeV, principally because of the 6=1/2 renormalon effect which were 
taken into account here and in Ref. Thus, the renormalon effect brings down the extracted central value of m,b 



7 We thank to M. Kalmykov for clarifications on this point. 
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TABLE VII: Rece ntly obtained values of (MS) mt, mass obtained from T sum rules or from spectrum of the T(1S) resonance. 
Wherever needed dl4ll35l ). the central mass values were adjusted to the common input central value a s (Mz) = 0.118. 



reference 


method 


order 


m b (m b ) (GeV) 


PP98 [431 


T sum rules 


NNLO 


4.21 ± 0.11 


MY98 [421 


T sum rules 


NNLO 


4.20 ± 0.10 


BS99 [67] 


T sum rules 


NNLO 


4.25 ±0.08 


H00 [57] 


T sum rules 


NNLO 


4.17 ±0.05 


KS01 [681 


T sum rules 


NNLO 


4.209 ± 0.050 


CH02 [69] 


T sum rules 


NNLO 


4.20 ± 0.09 


E02 [70] 


T sum rules 


NNLO 


4.24 ±0.10 


P01 [13J 


spectrum, T(1S) 


NNLO 


4.210 ±0.090 ±0.025 


BSV01 [55] 


spectrum, T(1S) 


NNLO 


4.190 ±0.020 ±0.025 


PS02 [35] 


spectrum, T(1S) 


N 3 LO 


4.349 ± 0.070 


L03 [14] 


spectrum, T(1S) 


N 3 LO 


4.19 ±0.04 


this work, Eq. H56H 


spectrum, T(1S) 


N 3 LO 


4.241 ±0.070 



TABLE VIII: Comparison of some of the toponium binding energy values E t i recently obtained in the literature. The first 
two values were correspondingly rescaled to our central input a 3 -value a s (Mz) = 0.1192, and rnt = 174.3 GeV. 

reference order E t i (GeV) 

PS02 [35] N 3 LO -3.065 ±0.157 (5=1,0) 

L03 [14] N 3 LO -3.21±0.15 (S=l) 

this work. Ea. l|I5afl N 3 LO -3.413 ± 0.153 (5=1) 

this work. Eg. ll49bl N 3 LO -3.481 ±0.163 (5 = 0) 



by about 150 MeV, but the separate evaluation/estimate of the ultrasoft contribution brings it up by about 50 MeV. 
The renormalon effect can also be understood from Fig. H^b) which suggests that (at \i « 3 GeV) the renormalon 
effect pushes upward the soft binding energy E b i[s) by about 300 MeV. We note that PS02 used pole mass mb in 
their N 3 LO TPS evaluation of the mass of the T(15) resonance before extracting the value of nib- 

Our results for the toponium binding energies are given in Eqs. I|46[) . I|48[) and i|49|) . in connection with Tables ITTT1 
and [El The result of Ref. [3, was E t{ « -3.08 ± 0.02 GeV (for S = 1), but the central value of a s used there 
was a s {M z ) = 0.1172. The result of Ref. H3 was E tI w -3.01 GeV, using the central value a s {M z ) = 0.1185. In 
Table IvTTTl we present our result (|49|) together with the results of these two references, in both cases rescaled to the 
common central a s value a s (Mz) = 0.1192. We see that in our case the toponium binding energies are significantly 
lower. This lowering is a combination of the renormalon b — 1/2 effect and of the ultrasoft effect. Fig.lHJa) suggests 
that the renormalon effect, in comparison to N 3 LO TPS, brings down the soft binding energy E t i{s) by about 150 
MeV (^ = 30 GeV) and 300 MeV (^ = 60 GeV). Further, the ultrasoft effect gBJ) brought down the binding energy 
by about 200 MeV. More specifically, when making our resummation with no separation of s and us part, and using 
the common renormalization scale fi = 50-60 GeV, would give results for the binding energy E t i higher by about 200 
MeV. The deviation of our result for E t i from the result of L03 in Table IVIIII can be explained principally with the 
ultrasoft effect, and the deviation from PS02 with combination of both the ultrasoft and the renormalon effect. We 
note that P02 used in their calculation of E t i the N 3 LO TPS with /i w 30 GeV and the pole mass m t . 

This lower binding energy E t i is then reflected in the value of the peak (resonance) position E res , - Eqs. (|55|l and 
Tables IVl and IVll For example, Ref. |3j| obtains for m t — 174.3 GeV [and central value a s {Mz) — 0.1192] the values 
E les . = 345.63 ± 0.16 GeV for 5 = 1 and 0, while we get the values 345.28 ± 0.26 GeV (5 = 1) and 345.21 ± 0.26 
GeV(5 = 0), i.e. lower by 350 and 420 MeV than [j^. In Ref. p3, NNLO results for E ICS . of several groups 
HH El E3- EI \IM were compared who used in their calculations various renormalon-free masses for the top quark. 
Their results were taken for the central input values a s {Mz) = 0.1190 and m t = 165.00 GeV, and are all around 
i?res. w 345.5 GeV, with variations, due to the renormalization scale ambiguity, being usually below 10 MeV. For 
these central input values of a s and m t , our results (J12J (see also Tables [V] and ED get modified to 347.34 ±0.26 GeV 
(5=1) and 347.27 ± 0.26 GeV (5 = 0), i.e., lower by about 200-300 MeV. 
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APPENDIX A: COEFFICIENTS FOR THE EXPANSION OF THE SOFT BINDING ENERGY 

We write down here the explicit coefficients fj of the expansion 126|) for the soft part of the ground state binding 
energy. The logarithms appearing in these expressions involve three scales {/i, jl,m q and 7l(/i) = (4/3)m 9 a s (pi)] 

^KS))' 



The coefficients fj are 



fi = ^(35 + 22L a - llLp - 11L 2 ) + |(-11 - 6Lj + 3L„ + 3L 2 )n f . (A2) 



H = f^+f^nf + f^n} , (A3a) 
/ 2 (0) = [381.674 + 90.75L 2 + 30.25^ + Li (246.417 - 121L M - 60.5L 2 ) - 48.5L 2 

+L^(-205.25 + 60.5L 2 ) - 11.69735(5+1)] , (A3b) 
/ 2 (1) = [ - 42.7469 - ILL? - 3.66667L^ + L M (26.6944 - 7.33333L 2 ) + 6.80556L 2 

+Li (-33.0556+ 14.6667L M + 7.33333L 2 )] , (A3c) 
/ 2 (2) = [1.16286 + (3/9)i? + (1/9)L* + L x (l - (4/9)L M - (2/9)L 2 ) 

+L M (-0.814815 + (2/9)i 2 ) - 0.185185i 2 ] . (A3d) 



h = /3 ( ° ) +/3 (1) n/+/fV / +/ 3 (3 M' (A4a) 
/ 3 (0) = [6726. 11 + 665. 5^- 166.375^ (40.8024 + (-10.5992 + L M )L M ) 

+L?(2381.5 - 1497.38L M - 499.125L 2 ) - 871.429L 2 - 499. 125(-1. 8843 + L M )L MJ L 2 
-201.438£| + Li (7457.17 - 497.292L 2 + L M (-4346.38 + 998.25L M + 998.25L 2 )) 
-257.341(0.211191 + L t - 0.75L M - 0.25L 2 )5(5 + 1) 

-61.4109 (-6.13937 + 5(5 + 1)) In (a a (^ a )) + 440.172 1h(k) + 2a 3 /4 3 ] , (A4b) 
= [ - 1274.33 - 1277.92Li - 471.125L? - 121Z? + 1182.32i |t + 843.6672,1^ + 272.25L?£ M 
-335.813L* - 181.5L X L^ + 30.25L 3 + 124.501L 2 + 108.361LiL 2 + 90.75L 2 L 2 
-186.708L M L 2 - 181.5iii^i 2 + 90.75-L^i 2 + 36.7292ZI 

+ (4.06858 + 15.5964Li - 11.6973^ - 3.8991L 2 )5(5 + 1)] , (A4c) 
/ 3 (2) = [70.8992 + 70.2453Li + 28.9722L? + 7.33333L? - 65.9925^ - 51.6667^1^ 
-16.5L?L M + 20.5972L 2 , + \\L x h\ - 1.83333L 3 - 5.19388L 2 - 6.57407LiL 2 

-b.bL\L 2 + 10.9167L JU L 2 + \\L X L^L % - 5.5^^ - 2.09722L 2 ] , (A4d) 
/f 5 = [ - 1.21475 - 1. 21714Li - (5/9)L? - 0.148148L? + 1.16286L M + LiL^ + (1/3) L^L^ 
-0.407407L 2 - (2/9)Lii^ + 0.037037L 3 + 0.0542857L 2 + (l/9)Ia£ 2 + (1/$)L\L 2 
-0.185185L M L 2 - (2/9)LiL M L 2 + (1/9)L 2 ! L 2 + 0.037037L^] . (A4e) 



22 



M. Peter, Phys. Rev. Lett. 78 (1997) 602 arXiv:hep-ph/9610209 ; Nucl. Phys. B 501 (1997) 471 arXiv:hep-ph/9702245 . 
Y. Schroder, Phys. Lett. B 447 (1999) 321 arXiv:hep-ph/9812205 ; Nucl. Phys. Proc. Suppl. 86 (2000) 525 
arXiv:hep-ph/9909520 . 

B. A. Kniehl, A. A. Penin, V. A. Smirnov, and M. Steinhauser, Nucl. Phys. B 635 (2002) 357 |arXiv:hep -ph/0203166|. 
N. Brambilla, A. Pineda, J. Soto and A. Vairo, Phys. Rev. D 60 (1999) 091502 arXiv:hep-ph/9903355 ; 
B. A. Kniehl and A. A. Penin, Nucl. Phys. B 563 (1999) 200 arXiv:hep-ph/9907489 . 
R. Tarrach, Nucl. Phys. B 183, 384 (1981). 

N. Gray, D. J. Broadhurst, W. Grafe and K . Schilcher, Z. Phys. C 48 (1990) 673; L. V. Avdeev and M. Y. Kalmykov, 
Nucl. Phys. B 502, 419 (1997) arXiv:hep-ph/9701308 . 

K. G. Chetyrkin and M. Steinhauser Nucl. Phys. B 573 (2000) 617 |arXiv:hep-ph/991 1434| ; K. Melnikov and T. v. Rit- 
bergen, Phys. Lett. B 482 (2000) 99 arXiv:hep-ph/9912391 . 

A. H. Hoang, M. C. Smith, T. Stelzer and S. Willenbrock, Phys. Rev. D 59 (1999) 114014 arXiv:hep-ph/9804227 . 
N. Brambilla, A. Pineda, J. Soto and A. Vairo, Nucl. Phys. B 566 (2000) 275 arXiv:hep-ph/9907240 . 
M. Beneke, Phys. Lett. B 434 (1998) 115 [arXiv:hep-ph/980424T1 . 

M. Neubert and C. T. Sachrajda, Nucl. Phys. B 438, 235 (1995) arXiv:hep-ph/9407394 ; M. E. Luke, A. V. Manohar 
and M. J. Savag e, Phys. Rev. D 51, 492 4 (1995) [arXiv:hep-ph /9407407]; A. G. Grozin and M. Neubert, Nucl. Phys. B 
508, 311 (1997) [arXiv:hep-ph /970731 8|; F. Campanario A. G. Grozin and T. Mannel, Nucl. Phys. B 663, 280 (2003) 
[Erratum-ibid. B 670, 331 (2003)J arXiv:hep-ph/0303052 . 
A. Pineda, JHEP 0106 (2001) 022 arXiv:hep-ph/0105008 . 
T. Lee, JHEP 0310, 044 (2003) |arXiv:hep- ph/0304185 . 

I. I. Bigi, M. A. Shifman, N. G. Uraltsev and A. 1. Vainshtein, Phys. Rev. D 50 (1994) 22 34 larXiv:hep-p h/9402360 . 

A. H. Hoang, Z. Ligeti and A. V. Ma nohar, Phys. Rev. Lett. 82, 277 (1999) |arXiv:hep-ph/9809423| ; Phys. Rev~D 59, 

074017 (1999) arXiv:hep-ph/9811239 . 

0. I. Yakovlev and S. Groote, Phys. Rev. D 63, 074012 (2001) arXiv:hep-ph/0008156 . 
T. Lee, Phys. Rev. D 67, 014020 (2003) arXiv:hep-ph/0210032 . 

1. Caprini and J. Fischer, Phys. Rev. D 60, 054014 (199 9) |arXiv:hep-ph/9811367| . 

G. Cvetic and T. Lee, Phys. Rev. D 64, 014030 (2001) arXiv:hep-ph/0101297 ; G. Cvetic, T. Lee and I. Schmidt, Phys. 
Lett. B 520, 222 (2001) arXiv:hep-ph/0107069 ; 

N. N. Khuri, Phys. Le tt. B 82, 83 (1979); G N. Lovett-Turner and C. J. Maxwell, Nucl. Phys. B 452, 188 (1995) 
arXiv:hep-ph/9505224 ; P. Ball, M. Beneke and V. M. Braun, Nucl. Phys. B 452, 563 (1995) arXiv:hep-ph/9502300 . 

G. Cvetic, Phys. Rev. D 67, 074022 (2003) |arXiv:hep-ph/0211226|. 

M. Beneke and V. M. Braun, Nucl. Phys. B 426 (1994) 301 arXiv:hep-ph/9402364 ; 
M. Beneke, Phys. Rept. 317 (1999) 1. 

M. Beneke, Phys. Lett. B 344, 341 (1995) |arXiv:hep -ph/9408380] . 
P. M. Stevenson, Phys. Rev. D 23, 2916 (1 981). 

G. Cvetic, Phys. Lett. B 486, 100 (2000) arXiv:hep-ph/0003123 ; G. Cvetic and R. Kogerler, Phys. Rev. D 63, 056013 
(2001) arXiv:hep-ph/0006098 . 

O. V. Tarasov, A. A. Vladimi rov and A. Y. Zharkov, Phys. Lett. B 93, 429 (1980); S. A. Larin and J. A. Vermaseren, 

Phys. Lett. B 303, 334 (1993) arXiv:hep-ph/9302208 . 

T. van Ritbergen, J. A. Vermaseren and S A. Larin, Phys. Lett. B 400, 379 (1997) |arXiv:hep-ph/9701390|. 

J. R. Ellis, I. Jack, D. R. Jones, M. Karliner and M. A. Samuel, Phys. Rev. D 57, 2665 (1998) arXiv:hep-ph/9710302 . 
V. Elias, T. G. Steele, F. Chishtie, R. Migneron and K. B. Sprague, Phys. Rev. D 58, 116007 (1998) arXiv:hep-ph/9806324 . 
T. Lee, Phys. Rev. D 56 (1997) 1091 arXiv:hep-th/9611010 ; T. Lee, Phys. Lett. B 462 (1999) 1 |arXiv:he p-ph/9908225 . 

G. Cvetic, arXiv:hep-ph/0309262 to appear in J. Phys. G. 

G. Cvetic, C. Dib, T. Lee and I. Schmidt, Phys. Rev. D 64, 093016 (2001) arXiv:hep-ph/0106024 . 

A. A. Penin and M. Steinhauser, Phys. Lett. B 538 (2002) 335 arXiv:hep-ph/0204290 . 

W. Fischler, Nucl. Phys. B 129 (1977) 157; A. Billoire, Phys. Lett. B 92, 343 (1980); for an application of a nonperturbative 
method see P. Gaete, Phys. Lett. B 515 (2001) 382 arXiv:hep-th/0107030 . 
T. Appelquist, M. Dine and I. J. Muzinich, Phys. Rev. D 17 (1978) 2074; 

B. A. Kniehl, A. A. Penin, M. Steinhauser and V. A. Smirnov, Phys. Rev. D 65, 091503 (2002) arXiv:hep-ph/0106135 . 
S. N. Gupta and S. F. Radford, Phys. Rev. D 24 (1981) 23 09; Phys. Rev. D 25 (19 82) 3430. 

S. Titard and F. J. Yndur ain, Phys. Rev. D 49 (1994) 6007 |arXiv:hep-p h/9310236 ; A. Pin eda and F. J. Y ndurain, Phys. 
Rev. D 58 (1998) 094022 arXiv:hep-ph/9711287 ; Phys. Rev. D 61 (2000) 077505 arXiv:hep-ph/9812371 . 
A. Czarnecki and K. Melnikov, Phys. Rev. Lett. 80, 2531 (1998) arXiv:hep-ph/9712222 ; Phys. Rev. D 65, 051501 (2002) 
arXiv:hep-ph/0108233 ; A. H. Hoang and T. Teubner, Phys. Rev. D 58, 114023 (1998) arXiv:hep-ph/9801397 ; J. H. Kiihn, 
A. A. Penin and A. A. Pivovarov, Nucl. Phys. B 534, 356 (1998) arXiv:hep-ph/9801356 . 
K. Melnikov and A. Yelkhovsky, Phys. Rev. D 59, 114009 (1999) arXiv:hep-ph/9805270 . 

A. A. Penin and A. A. Pivovarov, Phys. Lett. B 435, 413 (1998) |arXiv :hep-ph/9803363 ; Nucl. Phys. B 549, 217 (1999) 
arXiv:hep-ph/9807421 . 

F. A. Chishtie and V. Elias, Phys. Lett. B 521 (2001) 434 |arXiv:hep-ph/0107052|. 

M. Beneke and V. A. Smirnov, Nucl. Phys. B 522 (1998) 321 |arXiv:hep-ph/97H*39T] ; V. A. Smirnov, Phys. Lett. B 465, 



23 



226 (1999) arXiv:hep-ph/9907471 . 

Y. Kiyo and Y. Sumino, Phys. Lett. B 496, 83 (2000) arXiv:hep-ph/0007251 ; Phys. Rev. D 67, 071501 (R) (2003) 
arXiv:hep-ph/0211299 . 

S. Recksiegel and Y. Sumino, Phys. Lett. B 578, 369 (2004) arXiv:hep-ph/0305178 ; S. Recksiegel and Y. Sumino, Phys. 
Rev. D 67, 014004 (2003) arXiv:hep-ph/0207005 . 

G. Cvetic and R. Kogerler, Phys. Rev. D 66, 105009 (2002) arXiv:hep-ph/0207291 . 

K. Kajantie, M. Laine, K. Rummukainen and Y. Schroder, Phys. Rev. D 67, 105008 (2003) arXiv:hep-ph/0211321 ; 
J. P. Blaizot, E. Iancu and A. Rebhan, Phys. Rev. D 68, 025011 (2003) arXiv:hep-ph/0303045 . 

U. Aglietti and Z. Ligeti, Phys. Lett. B 364 (1995) 75 arXiv:hep-ph/9503209 . 

K. G. Chetyrkin, B. A. Kniehl and M. Steinhauser, Phys. Rev. Lett. 79, 2184 (1997) |arXiv:hep- ph/9706430]. 

G. Grunberg, Phys. Lett. 95B, 70 (1980), HOB, 501(E) (1982); 114B, 271 (1982); Phys. Rev. D 29, 2315 (1984); 

A. L. Kataev, N. V. Krasnikov, and A. A. Pivovarov, Nucl. Phys. B198, 508 (1982); A. Dhar and V. Gupta, Phys. Rev. 
D 29, 2822 (1984); V. Gupta, D. V. Shirkov, and O. V. Tarasov, Int. J. Mod. Phys. A 6, 3381 (1991). 

M. B. Voloshin, Nucl. Phys. B 154, 365 (1979); H. Leutwyler, Phys. Lett. B 98, 447 (1981); M. B. Voloshin, Sov. J. Nucl. 
Phys. 36, 143 (1982) [Yad. Fiz. 36, 247 (1982)]. 

B. L. Ioffe and K. N. Zyablyuk, Eur. Phys. J. C 27, 229 (2003) |arXiv:h ep-ph/0207183] 

N. Brambilla, Y. Sumino and A. Vairo, Phys. Rev. D 65, 034001 (2002) arXiv:hep-ph/0108084 . 
D. Eiras and J. Soto, Phys. Le tt. B 491, 101 (2000) |arXiv:hep -ph/0005066T~ 

A. H. Hoang, arXiv:hep-ph/0008102 

K. Hagiwara et al. [Particle Data Group Collaboration], Phys. Rev. D 66, 010001 (2002), and 2003 partial update for 
edition 2004 (URL: [http: / /pdg.lbl.gov). 

V. S. Fadin and V. A. Khoze, JETP Lett. 46, 525 (1987) [Pisma Zh. Eksp. Teor. Fiz. 46, 417 (1987)]. 
A. H. Hoang et al, Eur. Phys. J. directC 2, 3 (2000) |arXiv:hep-ph/0001286| . 

A. A. Penin and A. A. Pivovarov, Nucl. Phys. B 550, 375 (1999) |arXivT!ie p-ph/9810496 , order of NRQCD," Phys. Atom. 
Nucl. 64, 275 (2001) [Yad. Fiz. 64, 323 (2001)] arXiv:hep-ph/9904278 . 

M. Bohm, H. Spiesberger and W. Hollik, Fortsch. Phys. 34 (1986) 687; F. Jegerlehner, M. Y. Kalmykov and O. Veretin, 
Nucl. Phys. B 658, 49 (2003) |arXiv:hep-ph/0212319|. 

F. Jegerlehner and M. Y. Kalmykov, arXiv:hep-ph/0308216 

B. V. Geshkenbein, B. L. Ioffe and K. N. Zyablyuk, Phys. Rev. D 64, 093009 (2001) arXiv:hep-ph/0104048 . 
R. Barate et al. [ALEPH Collaboration], Eur. Phys. J. C 4, 409 (1998). 

S. Bethke, Nucl. Phys. Proc. Suppl. 121, 74 (2003) |arXiv:hep-ex/0211012|. 

M. Beneke and A. Signer, Phys. Lett. B 471, 233 (1999) |arXiv:hep-ph/9906475| . 

J. H. Kuhn and M. Steinhauser, Nucl. Phys. B 619, 588 (2001) [Erratum-ibid. B 640, 415 (2002)] arXiv:hep-ph/0109084 . 

G. Corcella and A. H. Hoang, Phys. Lett. B 554 , 133 (2003) |arXiv:hep-p h/0212297|. 
M. Eidemiiller, Phys. Rev. D 67, 113002 (2003) arXiv:hep-ph/0207237 . 

A. H. Hoang and T. Teubner, Phys. Rev. D 60, 114027 (1999) arXiv:hep-ph/9904468 . 
K. Melnikov and A. Yelkhovsky, Nucl. Phys. B 528, 59 (1998) arXiv:hep-ph/9802379 . 

O. I. Yakovlev, Phys. Lett. B 457, 170 (1999) arXiv:hep-ph/9808463 . 

T. Nagano, A. Ota and Y. Sumino, Phys. Rev. D 60, 114014 (1999) arXiv:hep-ph/9903498 . 



